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Abstract. The accuracy of the link of the proper motion system of astrometric satellite missions like AMEX and GAIA is 
discussed. Monte-Carlo methods were used to simulate catalogues of positions and proper motions of quasars and galaxies to 
test the link. The main conclusion is, that future satellite missions like GAIA may be "self-calibrated" by their measurements of 
QSOs, while additional measurements from radio stars or HST-data are needed to calibrate the less deep reaching astrometric 
satellite missions of AMEX type. 
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1. Introduction 

Since 1998 the reference frame for astrometric measure- 
ments is repre sented by the I nternational Celestial Reference 
Frame (TCRF. lMa et"all ll998). which is the realisation of the 
Intern ational Celestial Reference System (ICRS, lArias et all 
1 1995b . defined by the International Earth Rotation Service 
(IERS). The ICRF consists of 212 VLBI observed quasars as 
defining extragalactic radio sources. These objects are assumed 
to have no proper motion and no overall global rotation - the 
quasars represent a quasi inertial system. 

In the optical, the ICRS is r epresented today by the 
Hipparcos catalogue Jf.SAI 119971) . To be consistent, the 
Hipparcos catalogue had to be related to the ICRF. This means 
that it had to be assured that 1 st , the coordinate axes of the 
Hipparcos system are aligned with those of the ICRF and 2 nd , 
that the coordinate axes have no additional time dependent ro- 
tation component. Differences would show up as a global offset 
in positions and in a fictitious proper motion component. The 
transformation between the two systems is in general given by 
a rotation of the positional and proper motion system. This is 
the so-called link algorithm. 

For the Hipparcos satellite it was not possible to link the 
catalogue directly with extragalactic objects, because the only 
extragalactic object in the Hipparcos catalogue is the quasar 
3C 273 (HIP 6093 6). Consequently in termediate observations 
were needed JKovalevskv et all 1 19971) : VLBI, MERLIN and 
VLA observations of radio emitting stars, proper motions de- 
termined with the HST fine guidance system, and photographic 
observations with respect to extragalactic objects made fo r 
example in Bonn jGeffert et all 1 19971 iTucholke et all 1 19971) . 



The final link solutions were accurate to 0.6 mas for the 
posi tions and 0.25 ma s yr~' for the proper motion system 
dKovalevskv et all 1 9971) . 

The present study was initiated for the astrometric satellite 
mission AMEX 1 , now propo sed by the US Naval Observatory, 
and the ESA GAIA mission (Perrvman et al. 2001), with a fur- 
ther step in astrometric accuracy, foreseen to be launched no 
later than 2012. 

With the expected performance of AMEX some hundreds 
or even thousands of QSOs will be detectable, for GAIA even 
about 500 000 QSOs. If extragalactic objects are detectable it is 
possible to calibrate the proper motion system without further 
observations, so this would be a kind of self-calibration. 

The aim of this paper is to determine the accuracy of the 
link of the proper motions of the astrometric satellite missions. 
We present Monte-Carlo simulations of observations of QSOs 
with the performance expected for the AMEX satellite and also 
for the GAIA satellite mission. These simulations allow us to 
estimate the accuracy of a proper motion link using extragalac- 
tic objects. 

In Sect.|2]the link algorithm will be discussed. Some details 
about the implementation are described in Sect. 12.21 In Sect. [3] 
we present a search program for galaxies that could also be 
used for a proper motion link. In Sect. l4.ll we give a brief dis- 
cussion of the underlying model to simulate the observations 
and the testing is described in Sect. 14.21 Finally, in Sect. |5] we 
present the results of our simulations. 
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1 T he missions DIVA <Roser et alll997h and FAME <Horner et all 
are both cancelled, but the new satellite mission AMEX is a 
joint American-German merger of them. 
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2. The link algorithm 

In this section we will give a complete mathematical descrip- 
tion of the approximation of a small angle rotation and how it 
can be used to determine the components of the rotation vec- 
tors. 

The l ink algorithm has been used i n several publica- 
tio ns like fFroeschle & Kovalevskvl dl982l). lArias et al. ( 1988) 
or iLenhardt et alJ Jl997l) . In iLindeeren & Kovalevskv (1995) 
a theoretical study of the underlying exact model can be found. 
However, in most of these papers an approximation algo- 
rithm was used and the correction terms that will be described 
later (Eq. 1 1 31 were not taken into account, especially for the 
Hipparcos link program. We will show that with improved ac- 
curacy of astrometric data it is important to include these cor- 
rection terms. 

The basic idea of the algorithm is that in general the trans- 
formation between any two coordinate systems with the same 
origin is given by a rotation. If the coordinate axes of these 
systems are nearly coincident a first order approximation of the 
rotation can be used. It is required that the catalogues repre- 
senting these two systems have a high degree of internal con- 
sistency. This is the case for the ICRF, for the Hipparcos cata- 
logue, and it will be for future astrometric catalogues. 

2.1. The first-order approximation 

Let x, y and z be the axes of a cartesian coordinate system: 
x points towards the origin of right ascension, z towards th e 
celestial north pole and y — z x x. Following Murra vi lli 983l) . 
the coordinates of a unit vector r towards a celestial object in 
this coordinate system are given in general by 



! cos 6 • cos a 1 
cos 6 ■ sin a 
sin 5 



(1) 



The two orthogonal vectors p and q which point towards posi- 
tive right ascension and declination, respectively, are given by 



p = (zxr) 



- sin a 
cos a 




(2) 



Aa and A6 are observable quantities: differences between right 
ascension and declination of two astrometric catalogues. It then 
follows 



Ar = Aa cos 5 p + A5 q 

and with the Eqs. and Q, one easily finds 



(6) 



Aa cos 8 - Ar ■ p 

= - sin 5 cos a e\ - sin 6 sin a + cos 6 S3 (7) 



and 



AS — Ar ■ q — sin a e\ - cos a £2 



(8) 



where e, denotes the component of the vector s. For the proper 
motion components one finds 



A/j. a — —cos 6 cos a fis £1 + tan 5 sin or ji a s\ 

- cos 2 5 sin a us S2~ tan 5 cos a n a £2 

- tan 6 cos a u)\ - tan 6 sin a a>2 + W3 



and 

A/J.S — cos a ju„£i + sina /j a S2 
+ sin a oj\ — cos a a>2 



(9) 



(10) 



with a — n a , 5 - ids and e, = a>;. 

Finally the Eqs. @ - dlOt can be written in a matrix equa- 
tion that can be solved for the components of the rotation vector 
s and the spin vector 10 
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with the sub-matrices 



A,, = 



- sin 5 cos a - sin 6 sin a cos 6 
sin a - cos a 



(12) 



and 



and 



q = r x p = 



sin 6 ■ cos a ^ 
- sin 6 ■ sin a 
cos S 



(3) A 21 = 



where the brackets ( ) denotate vector normalisation. 

Any small rotation about an arbitrary axis e with value 
|e| <K 1 is given by 



Ar - exr 



(4) 



with the rotation vector s = \e\e. As shown in lMurravl lfl98.T) 
the vector e can be split into a rotation Aa about z and a rotation 
-A5 about p: 



e = Aa z - A6 p 



(5) 



sin 6 sin a fi a 
cos a n a 



cos 6 

sin 6 cos a ji a 
sin a /u a 





0) 



(13) 



The Equation (II 31 is important: this term corrects for the ori- 
entation difference of the vectors p and q on slightly different 
places on the celestial sphere. This is important because these 
are the directions in which the proper motion components are 
measured. 

For each celestial object one gets the 4 observational prop- 
erties and the matrix in Eq. (lilt which can be concatenated to 
a single matrix equation. For values AX we use the convention 
(reference - satellite) AX = X K f - X s . dt to calculate the matrix 
we use the reference catalogue. 
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2.2. Solving the equation 

The matrix equation dl It can be solved by a weighted linear 
least-squares method. We use the singular value decomposi- 
tion (SV D) because this alg orithm is numerically the most ap- 
propriate dPress et al.ll2002h . For the SVD, vector, and matrix 
operations we use the C++ library newmatlUA by R. Davies 2 . 
Details about the software implementation can be found in 
iMetTl ll2003ll . 

3. Galaxies for a proper motion link 

In principle galaxies are, besides quasars, also good candidates 
for a proper motion link if they are sufficiently far away and 
well detectable for the instrument. Ideal candidates would be 
galaxies that have a star-like core region for which the position 
can be measured by the satellite instrument. 

While - as mentioned above - only one extragalactic ob- 
ject can be found in the Hipparcos catalogue, the Tycho mea- 
surements should co ntain data for se veral galaxies. Since the 
Tycho-2 Catalogue JHe>g et alJ EoOO) has a completeness of 
90% at V = 1 1.5 mag, one would expect some bright galax- 
ies in this catalogue - however, not yet recognised as galaxies. 
On the other hand these objects may be important, since their 
positions from the Tycho-2 mission may be used as first epoch 
material. In this paragraph we have identified these galaxies 
and discuss the possibility of using them for the proper motion 
link of future satellite missions. 

3.1. Search program 

To identify galaxies in the Tycho-2 Catalogue we have com- 
pared the positions of all objects i n the Tycho-2 ma in catalogue 
and the Tycho-2 supplement 1 ( Hci g et all 120021) with cata- 
logues of galaxies and quasars available at the CDS FTP-Server 
( ftp://cdsarc .u-strasbg . fr/). The most important were 
the Updated Zwicky Catalogue (UZW), the Catalogue of 
principal galaxies (PGC), the Third Reference Catalogue of 
Galaxies (RC3) and the Arcsecond positions of UGC galaxies 
(AUGC) catalogue. 

All objects found were further reviewed: we checked 
each single object with Digitized Sky Survey (DSS) im- 
ages, whenever possible in different colors. In almost all 
cases miss-identifications due to stars in front of galax- 
ies could be eliminated. All other objects were fur- 
ther checked with the NASA/IPAC Extragalactic Database 
(NED, http://nedwww.ipac.caltech.edu/! and Simbad 
( http : //simbad . u- strasbg . fr/Simbad/ 1. 

Finally the galaxies were classified in two categories: 
"galaxies in the Tycho-2 Catalogue" and "uncertain galaxies 
in the Tycho-2 Catalogue". Objects were classified as uncer- 
tain in the case of: galaxies in dense star fields, a comment 
in the database that a star may be in front of the galaxy, or a 
clearly visible positional offset in the DSS images. Galaxies 
found in the Tycho-2 Supplement 1 were not classified as cer- 
tain/uncertain. 

2 Source code available for free from 
http : / /www . robertnz . net/nmlO . htm 



^ Galaxies 
"uncertain Galaxies 




Fig. 1. Galaxies in the Tycho-2 Catalogue plotted in galactic 
coordinates. The points in the background are the positions of 
all PGC galaxies. We have found 1 16 galaxies, 35 classified as 
uncertain and 30 galaxies in the Tycho-2 Supplement 1. The 
quasar 3C 273 was found (again) too. 

3.2. Galaxies in the Tycho-2 Catalogue 

A total number of 181 galaxies were found: 116 galaxies in 
the Tycho-2 Catalogue, 35 uncertain galaxies in the Tycho-2 
Catalogue and 30 galaxies in the Tycho-2 Supplement 1. The 
tables with all galaxies are available at CDS and the celestial 
distribution of the galaxies is shown in Fig.^ Beside the galax- 
ies, the quasar 3C 273 that was included in the Hipparcos cata- 
logue was found as well. 

All galaxies found have a PGC number but not all were 
found by the position cross correlation with the PGC due to 
the improper positional data in the PGC (and other galaxy cat- 
alogues as well). So there may be more galaxies in the Tycho-2 
Catalogue which could not be identified by the applied method 
without more accurate galaxy catalogues. 

However the galaxies found are hardly usable as first epoch 
material for a proper motion determination since they have 
mean position errors of cr a * = 117 + 22 mas and Ws 
1 27 ± 24 mas. With second epoch material e.g. from GAIA 
the resulting proper motion error would be of the order of 
cr^ ss 6 mas yr _1 . The large errors are mainly due to the fact 
that the galaxies are close to the limit of the Tycho-2 Catalogue 
with a mean magnitude of Vt = 12.5 mag. On the other side 
one would expect that a few months of observations with GAIA 
would produce results with a higher accuracy than it could be 
possible with help of the Tycho-2 Catalogue data. 

3.3. Astrometric Observations 

For three galaxies we made observations at our 1 m-Cassegrain 
telescope at th e Observa torium Hoher List with the CCD cam- 
era HoLiCam dSanner et al J 19981) in October 2002: PGC 7223, 
11404, and 65866. For these galaxies we determined the in- 
ter nal position accurac y with a software package developed 
by iGeffert et al.l II1997I) . by fitting a star profile to the galaxy 
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position. As reference we used the Tycho-2 Catalogue. For 
PGC 7223 and PGC 65866 the positions could be determined 
equally good as for stars in the fields with mean errors 5 
15 mas. Both galaxies have a very bright, star-like core region. 
For PGC 1 1404 we found a superimposed faint star but with a 
position clearly offset from the centre of the galaxy and from 
the Tycho-2 Catalogue position. 

The results give a hint that in principle the positions of most 
of the galaxies found could be determined quite well even if a 
star profile were fitted, as we did. 



4. Monte-Carlo Simulations: Models and tests 

4.1. Models for the simulations 

To analyse the possibility of linking the proper motion system 
of an astrometric satellite with QSOs we made Monte-Carlo 
simulations of observations of QSOs and for comparison also 
of observations of radio stars. The steps of the simulation were 
as follows: 
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A. Creation of a (synthetic) reference catalogue. 

B. Creation of a satellite catalogue by applying an exact rota- 
tion model. 

C. Applying an error model for the reference catalogue. 

D. Applying an error model for the satellite catalogue. 

E. Using the link algorithm to obtain the rotation and spin vec- 
tors. 

F. Comparison with the vectors used in step B. 

To create a reference catalogue we could either use a to- 
tally synthetic model that creates equally distributed apparent 
brightnesses and proper motions or make use of existing cata- 
logues to create lists of objects with adopted apparent bright- 
ness and proper motion distributions. The first model is useful 
for testing purposes, the second one for simulations of obser- 
vations. To apply a rotation to the synthetic reference catalogue 
we us ed the exact rotation equations JLindegren & Kovalevskvl 
1995). 

Errors could be simulated with the following models: no 
errors; normally distributed errors for testing; errors depend- 
ing on the apparent brightness of a simulated object should be 
used to simulate the performance of an astrometric satellite (see 
Sect. [5}; error distribution of existing catalogues could be used 
to simulate a realistic error distribution for example of radio ob- 
servations. The adopted errors were applied to the catalogues 
by assuming each single error to be normally distributed around 
zero. 3 

We allowed different parameters for the simulations: total 
number of simulations N s \ m , number of simulated catalogues 
N cat , and number of simulated objects iVobj- For each simula- 
tion we calculated the RMS of the difference between input and 
derived rotation and spin vector, <x(s,) and <r(u>i), which are a 
measure for the mean expected errors of each vector compo- 
nent. 

3 This does not imply that the total error distribution must be normal 
distributed! 



Fig. 2. Testing the algorithm performance: the three solid lines 
show the mean errors of the position link cr(e, ), the broken lines 
the mean errors of the proper motion link <x(a>,), when different 
parameters were varied: s and u> varied (dotted), only 10 varied 
(short dashed) and correction terms unused (long dashed lines). 
The errors of the position link cr(e,) are not affected by the 
variations of e or u>. 

4.2. Testing the algorithm 

We made detailed tests of the simulations and the algorithm. 
The accuracy of the link algorithm itself could be tested with 
simulations without errors. We varied the following parame- 
ters: 

- number of objects N ty 

- absolute value of the rotation and spin vectors |e| and |o>| 

- the absolute value of the simulated proper motion of objects 

Varying the number of objects, N h, does not lead to a change 
in accuracy. 

The effects of the correction terms of Eq. 1131 were tested 
by varying \s\ and \io\ with fixed mean proper motions, which is 
shown in Fig. [2 The error of the position link cr(e,) increases 
with increasing |e|. The value cr(<y,) increases similarly if |e] 
and \(o\ increase simultaneously. If in contrast the correction 
terms are not taken into account, the error <x(<y,) is larger by a 
factor x 500. The error cr{cjj) increases also if only the rotation 
vector |e] is varied whereas the value cr(e,) is not affected by 
any variation of \oj\. 

The different behaviour for small values of |e| and |w| is a 
numerical effect. Important for the numerical accuracy of the 
link algorithm are differences between two values, for exam- 
ple Aa. Whereas for the proper motions one has to store small 
differences of small values, for positions one has to store small 
differences of large values. The smallest absolute value |A<ar| 
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Fig. 3. Testing the error behaviour with normally distributed 
errors for the simulated reference and satellite catalogues. The 
mean proper motion error cr(w,) of the link depends on the size 
of the mean errors in the reference catalogue AR e f and the satel- 
lite catalogue Ac a t- The surface of error values is, of course, 
concave. 



that could be stored in a double precision floating point num- 
ber for a object with a > 15^3 is |Aa| = 2 3 " 52 w 4 • 1(T 7 mas. 4 

Also visible in the plot are little "bumps". These appear be- 
cause the RMS varies due to the Poisson-like distribution while 
a median would be more appropriate. However, in the plot we 
show the RMS because this it better comparable with the fur- 
ther results and shows the same basic behaviour as a median. 

Further tests showed that the mean errors are independent 
of the value of the mean proper motion /7 if the correction 
terms are used. If the simulated objects have no proper motion 
the results do not differ whether the correction terms are taken 
into account or not. 

We also made simulations with normally distributed errors 
for positions and proper motions of both catalogues. If we fix 
the mean errors of both catalogues and vary the number of link 
objects A^ bj the resulting error could, as one would expect, be 
fitted by a function oc 1/ JN bj . 

The relation between the resulting link mean error and the 
mean errors of the catalogue is shown in Fig. [5] We plotted the 
error cr(w,) depending on the size of the mean errors A of the 
simulated reference (ref) and satellite (cat) catalogues. The be- 
haviour is as follows: if we fix the mean error of one catalogue, 
for example A cat = mas yr , and increase the other one, the 
resulting error for the proper motion link cr(w,) increases lin- 
early. If we now fix A le f = 4 mas yr" 1 and increase A cat , one 
sees that the link error is further on dominated by the much 
larger error A cat until the value of A re f gets in the order of A cat . 
Finally, both contribute equally to the resulting error. 

If we assume normally distributed errors for both cata- 
logues a^ 2 -Test could be applied, but it fails if we use models 
for astrometric satellites (see next section), because the errors 
are no long er normally d istributed, which is a requirement for 
the ^ 2 -Test dBrandi 1999b . 

4 For a IEC/IEEE 64 bit floating point number, e.g. a double in C++ 
on i386 architecture. 




Wj [mas yr '] 



Fig. 4. Two components of the rotation vector to for a simu- 
lation of the proper motion link with quasars with a limiting 
magnitude of V max = 15 mag and the AMEX model. The com- 
ponents used (cross) are reproduced without bias but the stan- 
dard deviations are quite large with o-(a>i) = 0.64 mas yr -1 and 
cr(«2) = 0.58 mas yr" 1 . 

5. Simulations for astrometric satellite missions: 
AMEX like satellites and GAIA 

In lHirtd Jl998l) the relation between the apparent brightness 
of a star and the parallax accuracy for the DIVA satellite is 
presented, which then could be transformed in a position and 
proper motion accuracy. We fitted t wo exponential functions 
through the data points (Hirte 2002) to get a smooth relation 
and we use this data as estimation for AMEX like satel lites 
(AME X model). For GAIA we used data from lPerrvman et all 
and fitted a single exponential function (GAIA model). 
We then used these models to simulate a proper motion link 
with different kinds of objects. 

5.1. Simulations with quasars 

To simulate observations of quasars we used the cat- 
alog ue Quasars and Acti ve Galactic Nuclei ( QAGN) 
dVeron-Cettv & Veronl l200ll) to create the reference cata- 
logue with a luminosity distribution of quasars as realistic as 
possible. We choose different limiting magnitudes to deter- 
mine how the accuracy behaves with satellite performance. All 
proper motions and their errors were set to mas yr" 1 for the 
reference catalogue. This means that the resulting accuracy 
is only depending on the expected errors of the satellite 
observations. 

The parameters for the AMEX like satellite simulation 
were a limiting magnitude of V max = 15 mag (A^obj = 88 QSOs) 
and we made N slm = 10000 simulations. In Fig.|4]the results 
of one of those simulations are shown. It is clearly visible that 
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the values used to create the satellite catalogue (marked by the 
cross) are found on average but the results spread over a large 
range. 

The results of all simulations performed are given in Table 
^ The expected errors of the proper motion link for GAIA are, 
at the same limiting magnitude and the same number of objects, 
a factor of * 400 - 800 smaller than for the AMEX model, 
based on the expected GAIA performance llPerrvman et alJ 
120011) . 

The mean error of the proper motion link for the AMEX 
model is <x(<y ; ) = 0.61 mas yr" 1 at a limiting magnitude of 
V max = 15 mag. If we include more QSOs by using fainter 
objects the mean error decreases only slightly. This is due to 
the fact that simultaneously the error of the proper motion de- 
termination increases for fainter objects. So the profit of using 
more objects is essentially cancelled out by the worsened per- 
formance. For GAIA the error decreases steadily. 

To study the influence of the individual parameters we have 
also varied the number of QSOs and the length of the mission, 
both increased by a factor of 1.5. The results were as expected: 
if we increase the number of objects by a factor of 1 .5 the er- 
ror of the proper motion link decreases by a factor of Vl.5, 
whereas a 1 .5 times longer mission leads to a link error that is 
by a factor Vl.5 3 smaller jHirtell 998). 

In lKovalevskvH2003l) the effect of the curvature of the Solar 
system barycentric motion around the centre of our galaxy is 
described. This appears as an apparent proper motion of up to 
4 fias yr" 1 . Our simulations show that this affects the proper 
motion link by an error that is a factor of « 10% larger for all 
limiting magnitudes for the GAIA model if the data is not cor- 
rected for the aberration effect (which can be done before the 
link is applied). The apparent proper motion has no significant 
influence on the AMEX model with quasars. 

It has to be taken into account that the QAGN is not com- 
plete for faint objects. If we extrapolate the number of quasars 
brig hter than Vmxx — 1 6.5 mag in the SDSS Quasar Catalog I 
( Schneid er et al.ll2002l) we find that at this limiting magnitude 
about 57% ± 16% of the quasars are in the QAGN whereas only 
5.4% ±0.1% brighter than V max = 19 mag are in the QAGN. 
So for GAIA we underestimate the number of quasars in our 
simulations by a large factor. 

5.2. Simulations with galaxies 

As described in Sect. [3] we identified 116 galaxies in the 
Tycho-2 Catalogue. We used the position and magnitude data 
from the Tycho-2 Catalogue to simulate a proper motion link 
for AMEX like satellites with galaxies. All proper motions 
and their errors were set to mas yr" 1 for the reference cat- 
alogue, for the satellite errors we used the AMEX model. The 
resulting mean errors of the proper motion link were cr(w,) = 
0.16 mas yr" 1 . 

This result is a crude estimation. Firstly we made the as- 
sumption that the astrometric errors only depend on the mag- 
nitude but they may be affected by the source not being point- 
like. Secondly we used Tycho-2 magnitudes as estimation of 



Table 1. Expected mean errors <x(w,) of the proper motion link 
when QSOs were used for the simulations with the expected 
AMEX and GAIA performances. Proper motions and their er- 
rors were set to mas yr" 1 for the reference catalogues. 
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Table 2. Expected mean errors <x(w, ) of the proper motion link 
when radio-stars were used for the link with the AMEX model. 
We used the data of the 1 1 radio-stars observed with the VLBI 
and used for the Hipparcos link. In column 3 the errors of the 
VLBI observations are smaller by a factor of 2. 
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0.021 


0.014 


100 


0.020 


0.013 



the brightness of the core region. On the other side it is likely 
that more galaxies may be detected (see Sect. I3.2I >. 

5.3. Simulations with radio stars 

To compare the results of the simulations with quasars and 
galaxies with other methods w e made simulations with radio- 
emitting stars (Lestrade et alJI 19991) which were used for the 
Hipparcos link program. Radio stars are suitable link objects, 
because their proper motions may be determined with high 
accuracy within short time scales. For the construction of 
the synthetic referen ce catalogue the astrometric data from 
iLestrade et aD Jl999t) were used, magnitudes were taken from 
the Hipparcos catalogue. We varied the number of objects N bj 
and in a second run the VLBI errors were reduced by a factor 
of 2. 

The results are given in Table |2] The errors are mainly 
caused by the VLBI errors. The accuracy of the proper mo- 
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tion link could be increased by using more link objects which 
could be described by a function oc 1/ y/N as shown in Sect. 14. 21 
for normal distributed errors. 

6. Conclusions and outlook 

The main goal of the proper motion link of future astrometric 
satellite missions is to establish a proper motion system with 
an accuracy that is better than the accuracy of the single proper 
motions of one star. Our simulations show that for astrometric 
missions like GAIA this may be achieved by the use of QSOs 
alone. For astrometric satellite missions with AMEX like per- 
formance the use of quasars alone is not sufficient and galax- 
ies have to be included. However, an accurate proper motion 
link requires external data from VLBI observations of radio- 
emitting stars or HST observations of stars with respect to ex- 
tragalactic objects. 

Acknowledgements. We thank the referee Jean Kovalevsky for his 
useful comments as well as Ulrich Bastian and Klaas S. de Boer for 
useful discussions. 
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